Opportunities for future supernova studies of cosmic acceleration 
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We investigate the potential of a future supernova dataset, as might be obtained by the proposed 
SNAP satellite, to discriminate among different "dark energy" theories that describe an accelerating 
Universe. We find that 

many such models can be distinguished with a fit to the effective pressure-to-density ratio, w, of 
this energy. More models can be distinguished when the effective slope, dw/dz, of a changing w is 
also fit, but only if our 

knowledge of the current mass density, Sim , is improved. We investigate the use of "fitting functions" 
to interpret luminosity distance data from supernova searches, and argue in favor of a particular 
preferred method, which we use in our analysis. 



PACS Numbers : 98.80. Eq, 98.80. Cq, 97.60. Bw 

Increasing evidence that the Universe is accelerating 
confronts us with deep unresolved questions. As yet 
no compelling understanding of the acceleration has been 
achieved but many models have been proposed, typically 
introducing a form of "dark energy" to account for cos- 
mic acceleration. Further resolving the observational ev- 
idence for acceleration is certain to have a great impact 
on our understanding of the "dark energy" and of the 
underlying fundamental questions. 

Type la supernovae can be used as standard candles to 
infer the luminosity distance (di) as a function of redshift 
, and such data provide a key element in the case for 
cosmic acceleration. By analyzing a simulated dataset, 
as might be obtained by the proposed SNAP satellite [|j , 
we can test the ability of such experiments to distinguish 
among currently attractive models. Our work will be re- 
ported more extensively in [||. We note that where our 
calculations overlap we are in complete numerical agree- 
ment with the recent paper by Maor et al. [Q. However 
our perspective and emphasis is considerably different. 

We consider a future supernova dataset that has been 
converted to a table of effective magnitudes and redshifts 
of objects with a single fiducial absolute magnitude M. 
We will consider both statistical and systematic uncer- 
tainties in the magnitudes. Typically the redshifts are 
known to sufficiently high precision that their uncertain- 
ties can be ignored. 

The magnitudes and redshifts are related by the lumi- 
nosity distance according to 



m{z) = M 51og2?L, 



(1) 



where 



/■^ c luminosity 
V, Hod, = H„(l + z)j^^^ dz' (2) ^^(,) ^ ^ 



and M = M — 51ogi7o + 25. When nearby supernovae 
are used to determine AA the uncertainties in the Hubble 
constant drop out of Eqn. |l| It is through 'Dl{z) that 
specific models of the evolution of the universe enter the 



picture. If we had just a few well-specified models of dark 
energy to compare, then their predictions for Vl could 
be fit to the supernova dataset and their relative likeli- 
hoods determined. Even if a model had a free parameter, 
that parameter could be determined to some accuracy. 
However, it is not understood what fundamental physical 
principles specify the actual dark energy in the Universe, 
so we are left with a potentially infinite set of families of 
theoretical models to compare, each with a potentially 
large number of free parameters. No finite dataset will 
ever select among all of these theoretical models (see, for 
example, Maor et al.) 0|, who show nine models, from 
this infinite set, that cannot be distinguished by a pro- 
posed supernova dataset). 

In the face of this theoretical uncertainty it is useful to 
select a fitting function whose parameters can character- 
ize the m(z) dataset independent of any specific theory. 
By choosing a function that also fits most extant theories' 
Vl predictions well, with the fewest possible parameters, 
we can then use these parameters as indices that label the 
theories so they can be compared with the dataset. Al- 
though it is also desirable for the fitting parameters to 
represent physically meaningful concepts of the current 
meta-theories, this is not necessary. This approach has 
already been undertaken by several authors in the con- 
text of supernova datasets [^-|io|,0 . Each of these papers 
uses different types of functions to fit 'Dl{z). Here we 
describe a reason to choose one of these functions. 

Original reports of supernova data used (two- 
parameter) models based on adding a cosmological con- 
stant term to Einstein's equations. This approach was 
generalized in to an effective, but constant, equa- 
tion of state factor. Huterer and Turner modeled the 
luminosity distance according to the simple power law 
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Saini et al. ||T^ used the form|^ 



^ As we completed this work a preprint came out using a fit- 
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VLiz) = 2(1 



z — oi\J\ + z + a 
I3z + 7VI + z + 2 - a - 7 



(3) 



Our preferred fitting function is motivated by cos- 
mic acceleration driven by an extra energy component 
("quintessence" or "dark energy") with the pressure p 
and density p related by an equation of state p = wp. The 
cosmological constant case is reproduced when w = —1. 
For more general constant values of w this equation 
of state implies p cx (1 + z)'^'-™"'"^^ To allow for a z- 
dependence of w, it has been expanded in a power series 
in z (as in [^2[0]) or,equivalently, in (1 + z) 
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We assume a flat universe (theoretically preferred by in- 
flation and currently favored by CMB data). Thus the 
remaining parameter required to fix the cosmology is ilm, 
the density of ordinary (pressureless) matter today in 
units of the critical density. For these models we get 



Vl{z)=c{1 + z) I (1 + z')"'^' 



i7„i + r^Q (1 + z 
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exp 35:^ (l + ^')^-l 



dz' (5) 



where for a flat universe Qq = 1 — flm. Equation 5 takes 
on a simple form in terms of Wi, but throughout the rest 
of this letter we use the more intuitive quantity wt . 

In 1^ we survey essentially all the published 
quintessence models, and in every case the w-expansion 
gave a better flt than either Eqn. ^ or the power law 
expansion of d^, when taken to the same order of fit 
variables (see Fig. |l|). The quadratic order di-expansion 
(top dashed line) results in a worse fit than the linear 
order w-expansion (top solid line), and the quadratic or- 
der w-expansion (lower solid line) fits as well as the cubic 
order di-expansion (lower dashed line). (This result also 
holds when we weight the data points according to the 
dataset specification described below, although all the 
fits improve.) 

The w-expansion parameters also have an intuitive 
meaning in their own right. Given the uncertainty re- 
garding causes of cosmic acceleration, it might be rea- 
sonable to simply think of the w^'s as the fundamental 
parameters for now. They tell us something relatively 
straightforward about the nature of the dark energy. The 
same cannot be said of the a of M or the a, (3 and 7 of 



ting function with even more parameters ||ll| . We will discuss 
this function in Q 
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FIG. 1. Residuals from tliree fitting expansions to a given 
model (the periodic potential given in |l^]). The dashed lines 
are the quadratic and cubic polynomial di-expansion fit, the 
solid lines are the linear and quadratic w-expansion fit and 
the dotted line the fit according to Eqn. 

]lO| . We are fortunate that the twin criteria of intuitive 
interpretation of parameters and good fits to published 
models point to the same choice of fitting functions. It is 
worth noting however that we have not proven that the 
ui-expansion provides the best possible fits to current 
models. It is possible that an even better choice could be 
discovered in the future. 

Maor et al. also favor the same w-expansion formal- 
ism. However we do not agree completely with their dis- 
cussion of the alternative methods [p|jl^ . Despite claims 
to the contrary in all the approaches use "fitting func- 
tions" . If for example a new understanding of dark en- 
ergy emerged in which the Ci of |^ were fundamental 
parameters, then the w-expansion would provide a bad 
fit to those models and would be a bad choice for the fit- 
ting functions (whereas dL{z) = X]t=i '-i'^* would be an 
"exact expression"). 

So we have argued that the w-expansion offers a good 
means of interpreting supernova data. Now we apply this 
expansion to simulated future supernova datasets. 

For concreteness, we assume a future dataset similar to 
one proposed for the SNAP satellite project with SNe 
observed uniformly within four different redshift ranges 
with the following different sampling rates: In the first 
range from z — — 0.2 we assume that their are 50 ob- 
servations, in the second and largest redshift range from 
z = 0.2 - 1.2 there are 1800 SNe and in the two high 
redshift bins, z = 1.2 — 1.4 and z = 1.4— 1.7, there are 50 
SNe and 15 SNe observations respectively. The statistical 
error in magnitude is assumed to be CTmag = 0.15, includ- 
ing both measurement error and any residual intrinsic 
dispersion after calibration. It is worth noting that con- 
siderations such as those raised in could change this 
redshift-distribution strategy to further optimize the im- 
pact of the data. We construct Monte Carlo simulated 
datasets, first for a fiat, cosmological constant model with 
rim = 0.3 and fl^ = 0.7. Figure || shows an example of 
simulated data. The points show binned data with sta- 
tistical error bars, where we combine 56 SNe in each bin. 

Figure shows the magnitude difference. Am, between 
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several dark energy models (solid) and the Monte Carlo 
fiducial model, along with the simulated binned data for 
that model. To guide the eye we also show a grid of con- 
stant w models (dotted). Note how for the supergravity 
model (labeled "SUGRA"), adding the hnear term, 
wi , improves the fit considerably (as shown by the dashed 
line, mostly covered by the solid line), compared to the 
best constant fit with wq — 0.74. For the two exponen- 
tial model ("2EXP") the constant fit seems already 
sufficient. The line labeled "TOY" corresponds to a toy 
model with wq = —0.6 and wi = —0.8, which is not 
well represented by a constant fit. All the models have 

= 0.3. 
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FIG. 2. Simulated binned data, with solid curves from a su- 
pergravity inspired model (with the best A*' = 1 fit shown 
as a dashed line, mostly covered by the solid SUGRA line), a 
strongly evolving toy model with wo = —0.6 and wi = —0.8 
(TOY), and a model involving two exponentials Q {2EXP). 
A range of constant w models is shown as a grid of lines. The 
shaded region corresponds to cosmological constant models 
with 0.25 <Vlm< 0.35 and Oa = 1 - Sim- 

Figure |^ illustrates that these models can be easily 

differentiated from one another by the simulated data, 
and that for the toy model the contribution 

from wi looks large enough to be discriminated by the 
data. The targeted SNAP systematic uncertainties can 
be estimated by considering for example a linear drift up 
or down in the data (as a function of z) which spans one 
error bar height over < z < 1.5 (Note that a constant 
systematic uncertainty would not affect the results.) Fur- 
ther investigation is required to determine if these targets 
are realistic. When the targeted systematic uncertainties 
are taken into account the "2EXP" and w — ~1 models 
can no longer be discriminated. In we report more 
thoroughly on the models which can and cannot have 
wi differentiated from zero by such datasets. In most 
published models wi is below the level observable by su- 
pernovae alone. 

We shall see below that the effect of introducing wi 



in the fit for this dataset is to substantially increase 
the uncertainty in wq (as emphasized in [Q). This is 
because with the additional fitting parameter there are 
many more ways of fitting the simulated data to good ac- 
curacy. This fact has nothing to do with the ability of the 
(in this case simulated) dataset to discriminate among 
the models shown. The ability to discriminate among the 
models is a fundamental property of the dataset, whereas 
the choice of fitting function is a matter of choosing suit- 
able analysis techniques. As with any experimental data, 
for a given dataset each of the shown models is a member 
of a degenerate family of models that cannot be differ- 
entiated by the data. Opening up the space of fitting 
functions allows one to probe this degeneracy. Further- 
more, we show below how realistic improvements in our 
independent knowledge of can break this degeneracy 
to an interesting degree. 

Figure ^ shows joint likelihood contours in the f7„i — wq 
plane for the simulated data, when these are the only pa- 
rameters used to fit the data. We show also the shift due 
to a linearly drifting systematic error of ±0.02 mag per 
1.5 units in redshift. Clearly the data will discriminate 
among the three models shown. 
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FIG. 3. Separation of dark energy models (labeled as in 
Fig. 2) in the Vi^-wo plane. We show 39.3% (Ax^ = 1, 
shown as solid-line contours) and 68% (Ax^ = 2.3, shown as 
shaded contours) confidence regions which project on the axes 
to give marginalized Icr and 1.5o" errors respectively. The left 
and right flanking contours for each model show the effects of 
a drifting systematic error (see text). The w — —0.7 model 
has fim = 0.6, the other two models are exactly the same as 
their counterparts in Fig. 2. More examples are given in [^. 

We now analyze the same simulated data using a fit- 
ting function that includes wi. For this analysis (results 
shown in Table 1 and Figure ^), we consider different 
possible states of prior information about the value of 
f2m, ranging from relatively poorly known (crji^ ~ 0.15), 
roughly the case today, to well determined ((Tsi„ ~ 0.04), 
a possible goal within the next decade. We show in |^ 
that the systematic errors lead to roughly the same con- 
tribution as the statistical errors for the case of a per- 
fectly constrained = 0.0. We conclude from Fig. ^ 
and Table 1 that expanding the fitting function to include 
wi will become useful when the value of flm is better con- 
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prior crn,„ 


measurement (Jmag 






No Qrn prior; wi — 


0.15 


0.06 




0.15 


0.15 


0.15 


0.6 


0.05 


0.15 


0.06 


0.2 




0.09 


0.05 


0.12 


0.04 


0.15 


0.05 


0.16 


(fixed Qm) 


0.15 


0.03 


0.12 



TABLE I. Statistical measurement uncertainties on wo 
(i.e., Wtoday) and wi, given supernova magnitude measure- 
ment uncertainty, (Jmag , and a range of uncertainties, crn„ , in 
the independent prior knowledge of flrn- (As in Fig. 4, the 
supergravity model is used here as the example, but the other 
models give comparable results.) The systematic uncertain- 
ties are 0.01 in wo and 0.10 in wi. 



strained than it is today. This result adds to the case for 
producing an independent determination of f2„j. 
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FIG. 4. Separation of dark energy models (labeled as in 
Fig. 2) in the wq-wi plane. The shaded regions are 39.3% 
joint probabilities for the different Qrn priors and the dashed 
and solid lines the corresponding 68% probabihties. For the 
example of the supergravity (SUGRA) model, a range of in- 
creasingly larger contours represents the result using prior 
knowledge on fim with increasingly poorer uncertainty, an^ . 
The black dotted curve corresponds to the projected- Ict error 
contour obtained for more optimistic dataset specifications 
(o"mag ~ 0.09; see text) with a prior of 0.25 < fim < 0.35. For 
the example of the TOY model, fianking contours show the 
effects of a drifting systematic error (see text). 

It may be possible to reduce the statistical uncertain- 
ties even further, by increasing the sample size and/or by 
improving the intrinsic luminosity calibration for each su- 
pernova. Figure^ also shows the smaller likelihood region 
(dotted ellipse) which one gets when taking CTmag = 0.09 
and doubling the number of supernovae in each bin. 
This means (see Table 1) that the la statistical error 
on wi improves from a^i = 0.6 for the current con- 



straints on ftm (solid line in Fig. 4), to cr^j = 0.2 for 
future flm constraints (dashed line), to very tight con- 
straints of CTitij ~ 0.12 for the "improved" scenario, with 
twice as many more-tightly-calibrated supernovae and 
fmag — 0.09 (dotted line). 

Figure ^ is similar to Fig. 2 from j^] . (To compare ac- 
tual numbers for predicted uncertainties, note that Maor 
et al. quote the full-width range of 95%-confidence con- 
tours, i.e. four times a la uncertainty. Also, our differ- 
ent contour shapes are due to our choice of a Gaussian 
probability distribution for vs. their tophat distri- 
bution.) Maor et al. emphasize the uncertainties that 
emerge when the fitting function is expanded to include 
Wq and Wi given our current knowledge of flm- We dis- 
agree with this emphasis and focus instead on the great 
impact a SNAP-type experiment will have. Even if our 
determination of Qm were not to improve over the next 
decade, SNAP-class datasets would still be powerful dis- 
criminators among models. For example, such a dataset 
could easily discriminate between w = — 1 (a cosmolog- 
ical constant) as the source of cosmic acceleration and 
several interesting quintessence models (such as SUGRA 
or others discussed in |^). A result either way would 
have profound 

implications. Furthermore, significant improvements 
on the determination of Qm are expected over the next 
decade and we have shown how sufficient improvements 
would allow the extra fitting parameter wi to become 
an "observable" and help to further differentiate theories 

We have investigated different analysis methods for su- 
pernova datasets, and used our preferred method to ex- 
plore the prospects of further supernova searches. It is 
clear that a dataset of the sort proposed with SNAP ^ 
presents an exciting opportunity to constrain theories of 
cosmic acceleration. 
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